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This  paper  covers  numerical  simulation  results  for  an  in-stream  fuel  injector  in  a  Mach 
number  2.0  flow.  Three  fuel  injection  pylon  configurations  studied  are  a  basic  pylon,  a  ramp 
pylon,  and  an  alternating  wedge  pylon.  The  first  pylon  configuration  is  a  baseline.  The 
latter  two  configurations  introduce  streamwise  vorticity  into  the  flow  to  increase  mixing  ac¬ 
tion.  Computational  Fluid  Dynamic  simulations  with  ethylene  fuel  injection  are  presented 
and  evaluated.  Four  parameters  used  for  comparison  among  the  pylons  are  streamwise 
vorticity,  total  pressure  loss,  mixing  efficiency,  and  flammable  plume  extent.  An  overall 
energy  comparison  is  also  accomplished.  It  is  found  the  two  pylons  with  streamwise  vortic¬ 
ity  generation  are  better  mixing  devices  that  have  an  energy  advantage  over  the  baseline 
for  a  finite  distance  downstream  of  the  pylon. 

I.  Introduction 

Intrusive  fueling  techniques  for  scramjet  combustors  are  advantageous  in  distributing  fuel  uniformly 
over  a  combustor  cross-section.  The  performance  of  a  class  of  intrusive  in-stream  pylon  injectors  reported 
previously1  is  the  focus  of  this  paper.  The  current  work  provides  numerical  simulations  of  ethylene  fuel 
injection  from  three  geometric  variants  of  the  injector  concept  into  a  cold  airflow  at  a  Mach  number  of  2.0. 
The  three  pylon  variants  are  designed  to  introduce  varying  degrees  of  streamwise  vorticity  into  the  flow  to 
increase  fuel/ air  mixing.  The  pylons  are  compared  to  determine  which  geometry  results  in  the  best  fuel 
spreading  and  mixing  into  the  main  airflow.  Total  pressure  losses  are  also  of  concern,  so  the  spreading  and 
mixing  capabilities  of  each  pylon  are  considered  relative  to  the  pressure  loss  they  create. 

II.  Pylon  Configurations 

The  three  pylon  configurations  are:  the  basic  pylon  (Fig.  1),  the  ramp  pylon  (Fig.  2),  and  the  alternating 
wedge  pylon  (Fig.  3).  The  construction  of  each  pylon  incorporates  a  two-piece  forward  and  aft  area.  The 
forward  area  contains  a  plenum  common  to  all  pylon  configurations.  The  aft  area  contains  a  constant  angled 
compression  ramp  on  the  basic  pylon  and  a  hypermixer  geometry  on  the  other  two  pylons.  The  mating  of 
the  two  pieces  results  in  tall/thin  rectangular  injection  slots  created  by  the  difference  in  the  inner  width  of 
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the  forward  piece  and  the  outer  width  of  the  aft  piece.  The  common  parameters  for  all  pylons  are:  Height 
=  75  mm,  Length  =  103  mm,  Frontal  Blockage  Area  =  1215  mm2,  Fuel  Port  (Slot)  Area  =  57  mm2,  Front 
Wedge  Angle  =  14.7°,  and  Front  Wedge  Nose  Radius  =  1  mm. 

The  basic  pylon  configuration  embodies  the  two  fundamental  fuel  injection  strategies  of  1)  maximizing 
fuel-air  interface  area  by  using  a  rectangular  slot  injector  spanning  a  large  percentage  of  the  pylon  height, 
and  2)  injecting  fuel  upstream  of  the  pylon  aft  area  so  fuel  and  air  begin  mixing  prior  to  leaving  the  pylon 
surface.  Figure  1  shows  the  basic  pylon  on  a  pedestal.  The  aft  area  of  the  basic  pylon  is  a  compression  ramp 
at  10.6°  (half  angle)  to  the  main  airflow.  The  fueling  slots  operate  at  sonic  conditions.  A  plenum  in  the 
forward  area  of  the  pylon  provides  fuel  to  the  slots.  This  plenum  is  fed  by  an  opening  at  the  bottom  of  the 
pylon  connected  to  the  pedestal.  Fuel  injection  is  accomplished  from  a  backward  facing  step  that  protects 
the  injectant  from  the  main  flow  for  a  short  distance  prior  to  joining  the  main  airflow  around  the  pylon. 


Isometric  View 


Isometric  View 


The  ramp  pylon  shown  in  Fig.  2  includes  eight  compression  ramps  on  the  aft  area  of  the  pylon.  These 
compression  ramps  are  14.4°  to  the  main  airflow  with  8.3°  of  ramp  taper  (sweep).  Past  studies  have  shown 
that  sweep  on  a  hypermixer  compression  ramp  increases  streamwise  vortices  and  fuel/air  mixing  compared 
with  no  sweep. 2-4  The  compression  and  sweep  angles  were  chosen  to  produce  the  same  frontal  blockage  area 
as  the  basic  pylon  configuration. 

The  alternating  wedge  configuration  shown  in  Fig.  3  includes  eight  alternating  wedges  on  the  aft  area  of 
the  pylon.  Past  studies  have  indicated  alternating  wedge  geometries  also  produce  streamwise  vortices  and 
increase  fuel/air  mixing.5,6  The  wedge  geometries  have  a  23.6°  angle  and  attach  to  the  aft  area  of  the  pylon, 
itself  having  a  14.4°  compression  angle  to  the  main  airflow.  This  configuration  is  sized  to  produce  the  same 
frontal  blockage  area  as  the  other  two  pylons. 
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III.  Grid  Construction 


Grid  generation  is  accomplished  using  a  combination  of  GRIDGEN  and  SOLIDMESH.  Boundary  domains 
are  constructed  in  GRIDGEN  using  imported  surface  geometry  from  the  computer  aided  solid  body  design 
software  SOLIDWORKS.  These  meshed  surface  domains  are  imported  from  GRIDGEN  into  SOLIDMESH 
and  cells  are  constructed.  The  three  dimensional  grids  for  each  pylon  are  an  unstructured  hybrid  of  rect¬ 
angular/triangular  surface  domains  and  tetrahedral/pentahedral  cell  volumes.  The  viscous  tailoring  of  the 
grid  near  the  pylon  boundary  surfaces  is  accomplished  to  the  refinement  required  to  utilize  the  boundary 
layer  wall  functions  in  FLUENT  (30  <  y+  <  200).  The  initial  wall  spacing  and  grid  refinement  near  the 
pylon  boundary  surfaces  is  the  same  for  all  three  pylons.  The  exterior  walls  of  the  computational  domain 
are  inviscid  for  all  three  pylons  in  order  to  limit  grid  size,  hence  the  grid  is  not  refined  near  the  outer  walls. 
The  wake  region  of  all  three  pylons  is  populated  with  a  dense  cell  mesh  extending  290  mm  downstream  of 
the  pylons.  The  cell  volumes  in  the  wake  region  are  approximately  1  mm3.  The  grids  for  the  three  pylons 
vary  between  5.5  -  6.7  million  cells.  Figure  4  shows  a  typical  grid. 


Figure  4.  CFD  internal  grid. 


Only  one  grid  is  constructed  for  each  pylon.  The  grids  are  constructed  similarly,  with  the  most  attention 
and  refinement  given  to  the  wake  region  of  the  pylons.  A  grid  convergence  study  is  not  accomplished  to 
optimize  the  grids.  Any  comparison  of  these  results  to  experimental  data  needs  to  take  into  account  the 
inviscid  exterior  wall  of  these  grid  domains.  Wind  tunnel  walls  exhibit  viscous  effects.  CFD  solutions  will  be 
shown  similar  to  experimental  data  at  one  plane  location  downstream  of  the  pylons  in  a  subsequent  section. 

IV.  Flow  Solution 

FLUENT  is  the  Computational  Fluid  Dynamics  (CFD)  code  employed.  The  K -u>  turbulence  model  with 
SST  (shear  stress  transport)  is  used  with  a  multi-species  formulation  of  the  compressible/ viscous  Navier 
Stokes  flow  equations.  The  result  sought  is  a  second  order  steady  state  or  time  accurate  solution.  No 
combustion  is  modeled,  only  species  transport.  The  cold  airflow  inlet  and  ethylene  injection  conditions  are 
shown  in  Table  1.  These  conditions  result  in  a  unity  momentum  flux  ratio,  q  (Eq.  1). 
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_  {lPM'2)  ethylene  _  (^ethylene 
q~  (7  PM2)air  (pM2)air 


Table  1.  Flow  conditions  for  CFD  simulation. 


Property/Condition 

Cold  Airflow  Inlet 

Ethylene  Injection 

Ratio  of  Specific  Heats 

1.4 

1.25 

Molecular  Weight 

28.97 

28.05 

Mach  Number 

2.0 

1.0 

Velocity 

518  m/s 

314  m/s 

Mass  Flow 

10.78  kg/s 

0.04  kg/s 

Total  Temperature 

300  K 

300  K 

Static  Temperature 

167  K 

267  K 

Total  Pressure 

310  KPa 

320  KPa 

Static  Pressure 

39.7  KPa 

178  KPa 

Flow  Area 

25,161  miri2 

57  mm2 

Solution  convergence  is  determined  three  ways  for  a  steady  state  solution:  1)  the  drag  on  the  pylon  is 
stable,  2)  the  residuals  of  mass,  velocity,  and  energy  are  reduced  at  least  three  orders  of  magnitude  from 
the  original  values,  and  3)  the  spatially  integrated  total  mass  flow  and  ethylene  mass  flow  are  constant 
downstream  of  the  pylons.  Two  of  the  three  pylon  solutions  are  steady  enough  to  meet  all  three  criteria 
given  enough  iterations  of  localized  time-stepping.  The  third  pylon  exhibits  a  highly  unsteady  flow,  and  does 
not  meet  the  mass  flow  requirement  using  localized  time-stepping.  Global  time-stepping  is  implemented  to 
find  a  time  accurate  flow  solution  for  the  third  pylon. 
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V.  Data  Collection 


The  equivalent  diameter,  de,  is  the  diameter  of  the  rectangular  slot  fuel  port  area  if  it  were  to  form  a 
circle.  The  equivalent  diameter  for  the  fuel  port  area  (57  mm2)  is  8.52  mm.  The  axial  distance  downstream 
of  the  injection  slots  is  normalized  by  de.  The  dense  cell  wake  region  extends  34  de  downstream  of  the  pylon 
base  (aft  face),  and  the  base  of  each  pylon  is  40  mm  (4.7  de)  aft  of  the  injection  slots. 

Fourteen  downstream  positions  are  selected  for  reduction  of  five  data  parameters:  streamwise  vorticity, 
total  pressure  loss,  mixing  effectiveness,  flammable  plume  extent,  and  potential  work  change  ratio.  These 
eight  positions  correspond  to  5.9  de,  7.4  de,  8.9  de,  10.6  de,  12.1  de,  13.6  de,  15.3  de,  16.8  de,  18.3  de,  20 
de,  24.7  de,  29.4  de,  34  de,  and  38.7  de  aft  of  the  injection  slots.  Figure  5  shows  the  fourteen  downstream 
positions  chosen  for  data  reduction. 

Streamwise  vorticity,  u>x,  is  defined  as  the  axial,  or  x-component,  of  the  curl  of  the  velocity  vector.  The 
absolute  value  of  the  quantity  is  used  for  the  integration.  The  mass  averaged  integral  of  the  absolute  value  of 
axial  vorticity  is  the  measure  of  streamwise  vorticity,  f lx  (Eq.  2).  Instead  of  presenting  absolute  numbers  of 
streamwise  vorticity,  relative  magnitudes  are  computed,  normalized  by  the  maximum  vorticity  value  present 
in  the  data. 


I  ^ X  I 
\lx 


l(vxf,)» 

J  \ojx\  pud  A 
f  pudA 


(2) 


Total  pressure  loss,  A,  is  calculated  from  the  mass  averaged  integral  of  total  pressure.  The  mass  averaged 
integral,  divided  by  the  wind  tunnel  total  pressure  and  subtracted  from  unity,  is  the  pressure  loss  (Eq.  3). 
Total  pressure  losses  in  this  specific  problem  are  relatively  small  (all  three  pylons  block  only  4.8%  of  the 
axial  flow  path  area) . 


PT 


A 


/  PrpudA 
f  pudA 


1  - 


Pt 

~Pr 

J-  tunnel 


(3) 


Mixing  effectiveness  is  quantified  from  a  comparison  of  the  local  mass  fraction  of  fuel/air  to  that  for 
stoichiometric  combustion.  This  method  requires  the  injected  fuel  have  a  known  stoichiometric  fuel/ air 
ratio,  fs.  The  mass  fraction  of  fuel  in  air,  a,  can  be  calculated  from  fs  and  the  equivalence  ratio,  </),  with 
Eq.  4.  For  example,  with  fs  =  0.068,  the  mass  fraction  of  ethylene  in  air  at  stoichiometric  conditions  (</>  = 
1)  calculated  from  Eq.  4  is  as  =  0.064. 


a 

a 


mass  of  fuel 


mass  of  fuel  +  air 

4>fs 

4>fs  + 1 


(4) 


The  mixing  effectiveness  formulation  assumes  the  least  available  reactant,  fuel  or  air,  is  the  reactant 
considered  for  mixing.  The  other  reactant  is  considered  completely  mixed  already.  For  instance,  in  a  globally 
fuel-lean  situation  (which  is  the  case  here) ,  the  fuel  is  the  mixing  concern,  so  the  air  is  considered  completely 
mixed.  Equation  5  shows  the  mixing  effectiveness  quantity,  rym,  in  a  globally  fuel  lean  situation.7,8 


Vm 


difuel, mixed  _  f Aa=0  aR PudA 
di-fuel, total  -I.A,.,  r,  aPudA 

la,  a  <  as 

1  as(l  —  a)/(l  —  as),  a  >  as 


(5) 
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The  area  of  integration  extends  to  the  edges  of  the  fuel  plume  where  the  fuel  mass  fraction  drops  to 
zero,  or  practically,  to  some  minimum  detection  mass  fraction,  a  >  amin  (Omin  =  0.001).  If  the  local 
fuel  mass  fraction  is  less  than  the  mass  fraction  required  for  stoichiometric  combustion,  the  fuel  is  fully 
mixed.  If  the  local  fuel  mass  fraction  is  greater  than  the  stoichiometric  mass  fraction,  the  fuel  is  unmixed 
to  some  degree.  A  weighting  function,  seen  in  the  bottom  line  of  Eq.  5,  is  used  to  quantify  how  unmixed 
is  the  local  mass  fraction.  The  fully  mixed  numerical  value  is  one,  and  the  completely  unmixed  numerical 
value  is  zero.  This  mixing  effectiveness  does  not  complete  the  fuel  mixing  picture.  Mass  fractions  below 
stoichiometric,  considered  perfectly  mixed,  include  those  below  the  threshold  for  combustion.  In  addition, 
flammable  mixtures  exist  above  stoichiometric,  and  the  formulation  also  does  not  take  this  into  account.  For 
a  more  complete  picture,  one  must  look  at  the  flammable  plume  extent. 

The  flammable  plume  extent  is  observed  in  two  quantities.  First  is  the  fuel  plume  area  with  mass 
fractions  falling  within  the  flammable  region.  This  quantity  is  normalized  by  the  injection  slot  area.  Second 
is  the  fraction  of  fuel  plume  falling  within  the  flammable  region.  The  first  quantity  compares  the  overall 
flammable  plume  extent  between  pylons.  The  second  quantity  determines  the  most  efficient  place  to  ignite 
the  flammable  mixture,  at  the  point  where  the  fuel  plume  exhibits  the  largest  fraction  of  flammable  area. 
The  equivalence  ratio  flammability  limits  for  ethylene  are  used  (0.41  <  <j>  <  5. 5). 9  Consequently,  from  Eq. 
4,  the  range  of  flammable  mass  fractions  of  ethylene  fuel  in  air  is  0.027  <  a  <  0.272.  The  flammable  plume 
area,  FPa,  and  the  flammable  plume  fraction,  FPf,  are  shown  in  Eq.  6.  The  fuel  injector  area,  Ae,  is  5.7E-5 
m2,  and  the  fuel  plume  area,  Ap,  is  the  area  of  fuel  having  mass  fractions  greater  than  0.001. 


FPa 


FPf 


Af 

Ae 

Af 


(6) 


Figure  6.  Fuel  plume  visual. 


A  measure  is  required  to  assess  the  overall  performance  of  the  pylons.  There  is  energy  loss  due  to  reduction 
in  flow  total  pressure.  There  is  energy  gain  from  heat  release  due  to  fuel  mixing.  The  combination  of  this 
loss  and  gain  is  an  energy  change  we  shall  call  the  work  potential  change  of  the  combustor  flow.  Energy  loss 
in  a  scramjet  combustor  normally  includes  drag  losses  (pressure  and  viscous)  and  heat  release  losses.  During 
heat  release  in  a  scramjet  combustor,  total  pressure  and  hence  Mach  number  are  reduced,  decreasing  the 
work  potential  of  the  combustor  flow.  Drag  losses  also  appear  as  a  total  pressure  loss,  lowering  the  Mach 
number  and  the  work  potential  of  the  combustor  flow.  Total  pressure  losses  are  entropy  gains,  as  derived 
from  the  Gibbs  equation  and  first  law  of  thermodynamics  (Eq.  7).  Entropy  gain  results  in  a  reduction  in 
energy  and  work  potential  of  the  combustor  flow. 
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Energy  gain  from  heat  release  in  the  combustor  can  be  expressed  through  the  first  law  of  thermodynamics 
in  Eq.  8.  Heat  addition  is  an  augmentation  in  work  potential  of  the  combustor  flow.  The  mixing  efficiency 
parameter  multiplied  by  the  flammable  plume  fraction,  ?7mFPf,  is  substituted  for  the  combustion  efficiency 
that  should  be  present  in  the  equation.  Although  not  the  combustion  efficiency,  this  multiplied  quantity 
allows  for  some  calculation  of  the  heat  release,  rhf  is  the  mass  flow  of  fuel,  hf  is  heat  release  per  kilogram 
of  fuel.  This  equation  offers  some  measure  of  the  energy  increase  from  fuel  supplied  by  the  pylon  given  the 
data  available. 


AWq  =  A  Q  =  rhf7?mFPf/if  (8) 

The  summation  of  energy  gain  due  to  heat  addition  and  the  energy  loss  due  to  entropy  is  a  measure  of 
the  work  potential  change  of  the  combustor  flow  due  to  pylon  fueling  (Eq.  9).  This  method  omits  effects 
due  to  heat  release  such  as  basic  flow  field  changes,  large  pressure  drops,  and  vorticity/mixing  rate  changes. 
However,  in  the  absence  of  a  defined  combustor  geometry,  we  believe  this  approach  is  a  reasonable  alternative 
to  a  full  examination  of  a  scramjet  flow  path  incorporating  these  pylons. 


AW  =  AWl  +  AWq  =  mf  ??mFPf/if  -  mT^Rln  (9) 

The  work  potential  change  is  a  representative  quantity,  meant  only  to  compare  the  performance  of  each 
pylon,  not  to  predict  the  actual  work  potential  in  a  combustor.  In  presentation,  the  work  potential  change 
of  the  alternating  wedge  and  ramp  pylons  are  divided  by  the  work  potential  change  of  the  basic  pylon  to 
give  a  work  potential  change  ratio. 

VI.  Numerical  Simulation  Results 

The  basic  pylon  exhibits  highly  unsteady  flow  in  the  form  of  vortex  shedding.  Localized  time  stepping 
results  in  reduced  residuals  about  three  orders  of  magnitude  from  original  values,  and  a  steady  drag  on 
the  pylon,  but  the  total  mass  flow  downstream  of  the  pylon  is  not  constant.  The  total  mass  flow  decreases 
and  the  ethylene  mass  flow  increases  in  magnitude  in  the  downstream  direction.  This  is  a  violation  of  mass 
conservation.  A  time  accurate  solution  is  needed  for  this  pylon.  An  estimate  of  the  vortex  shedding  frequency 
is  accomplished  to  calculate  a  time  step  required  to  resolve  the  unsteadiness  in  the  wake  flow.  A  Strouhal 
number  of  0.25  is  assumed  as  a  typical  non-dimensional  shedding  frequency  downstream  of  bluff  bodies  in 
high  Reynolds  number  flows.  The  flow  velocity  is  518  m/s  and  the  width  of  the  pylon  base  is  0.018  m.  Given 
Eq.  10,  the  vortex  shedding  frequency  estimated  is  about  7200  Hz. 

st=SY  (10) 

A  time  step  of  one  microsecond  is  selected  to  resolve  this  shedding  frequency.  Monitoring  the  Mach 
number  of  three  points  in  the  wake  region  of  the  pylon,  the  flow  solution  is  advanced  in  time  until  a 
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consistent  periodic  pattern  emerges.  Figure  7  shows  the  Mach  number  over  time  at  one  of  the  three  points. 
The  same  periodic  frequency  is  observed  at  all  three  points,  approximately  6900  Hz  (St  =  0.24). 


Figure  7.  Mach  number  at  point  in  wake  region  over  time. 


The  solution  at  the  final  time  step  is  kept  for  analysis.  The  solution  is  then  advanced  one  half  period 
further,  and  that  solution  is  also  kept  for  analysis.  Data  collection  is  accomplished  at  the  eight  downstream 
planes  in  both  instantaneous  snapshots,  one  half  period  separated,  and  averaged.  As  it  turns  out,  all  the 
integrated  quantities  are  virtually  identical  in  the  two  time  snapshots. 

A  planar  cut  in  the  y-direction  of  the  Mach  number  half  way  up  the  basic  pylon  height  and  eight  planar 
cuts  in  the  x-direction  of  ethylene  mass  fraction  downstream  of  the  pylon  are  shown  in  Fig.  8.  Only  the  first 
of  the  time  solution  points  is  displayed.  The  other  time  point  is  very  similar,  where  the  vortex  shedding  is 
at  a  different  point  in  its  phase. 

The  ramp  and  alternating  wedge  pylons  both  exhibit  a  constant  mass  flow  in  their  wake  region  and  meet 
the  other  steady  state  solution  convergence  criteria  using  localized  time  stepping.  This  is  not  to  say  the 
flow  downstream  of  these  pylons  is  perfectly  steady,  just  steady  enough  to  satisfy  the  convergence  criteria 
selected.  It  appears  the  streamwise  vorticity  production  of  the  ramp  and  alternating  wedge  pylons  promotes 
a  more  steady  flow  in  the  wake  region.  Vortex  shedding  is  not  present  in  these  solutions.  The  same  planar 
cuts  of  Mach  number  and  ethylene  mass  fraction  are  displayed  for  the  ramp  and  alternating  wedge  pylons 
in  Figs.  9  and  10. 
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(a)  Mach  number. 


(b)  Ethylene  mass  fraction. 


Figure  8.  Basic  pylon. 
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VII.  Pylon  Comparisons 


In  general,  the  most  interesting  comparison  area  between  pylon  wake  flows  is  from  the  pylon  base  to 
approximately  20  de.  The  mass  flows  in  the  wake  region  of  the  three  pylons  are  compared  in  Fig.  11.  For 
the  two  steady  solutions,  the  total  mass  flow  and  ethylene  mass  flow  are  constant  in  the  wake  region.  For  the 
time  accurate  solution,  the  total  mass  flow  and  ethylene  mass  flow  in  the  wake  region  vary  slightly,  but  do 
not  on  average  increase  or  decrease  with  downstream  distance.  The  mass  flows  oscillate  around  the  correct 
values  of  10.82  kg/s  and  0.04  kg/s  respectively. 


Figure  12.  Vorticity  and  pressure  loss  comparison. 


The  streamwise  vorticity  and  pressure  loss  produced  by  the  three  pylons  are  compared  in  Fig.  12.  It  is 
difficult  to  make  direct  downstream  position  comparisons  between  the  basic  pylon  and  the  other  two  with 
the  slight  mass  flow  variations  in  the  basic  pylon  time  accurate  solution.  However  there  are  general  trends. 
The  general  trend  for  streamwise  vorticity  is  to  decrease  with  downstream  distance.  The  general  trend  for 
pressure  loss  is  to  increase  with  downstream  distance.  The  vorticity  present  immediately  downstream  of  the 
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ramp  and  alternating  wedge  pylons  is  increased  over  the  basic  pylon  due  to  the  vortical  inducing  geometries 
at  the  aft  end  of  those  pylons.  However,  the  increased  vorticity  of  these  geometries  is  observed  for  a  short 
distance.  In  this  narrow  window  of  space  is  the  opportunity  for  the  ramp  and  alternating  wedge  pylons  to 
mix  fuel  and  air  more  rapidly  than  the  basic  pylon.  Looking  past  the  near-held,  the  streamwise  vorticity 
downstream  of  the  basic  pylon  exceeds  that  of  ramp  pylon  and  the  alternating  wedge  pylon  as  downstream 
distance  increases  into  the  far-held.  The  pressure  loss  of  the  basic  pylon  exceeds  that  of  the  ramp  pylon 
as  downstream  distance  increases  into  the  far-held,  but  does  not  exceed  the  pressure  loss  of  the  alternating 
wedge  pylon. 

Overall,  the  increased  streamwise  vorticity  of  the  ramp  and  alternating  wedge  pylons  in  the  near-held 
wake  region  comes  with  added  pressure  loss  over  the  basic  pylon.  As  downstream  distance  increases  into  the 
far-held,  the  basic  pylon  comes  close  to  or  exceeds  the  ramp  and  alternating  wedge  pylons  in  vorticity  and 
pressure  loss.  The  alternating  wedge  pylon  has  larger  vorticity  and  pressure  loss  than  the  ramp  throughout 
the  wake  region. 

The  non-dimensional  drag  from  surface  pressures  is  the  drag  coefficient,  which  is  defined  in  Eq.  11.  The 
dynamic  pressure  multiplied  by  the  pylon  frontal  area  is  134.9  N  in  this  cold  how  for  all  pylons  since  each 
has  the  same  frontal  area  by  design.  The  drag  coefficients  are  presented  in  Table  2. 


Cd  = 


Drag 

-j  pu2(Area) 


Drag 

7^  PM2  (Area) 


(11) 


Table  2. 

Pylon 

drag. 

Pylon 

cd 

Drag  (Newtons) 

Basic 

0.61 

82 

Ramp 

0.59 

80 

Alternating  Wedge 

0.66 

89 

The  mixing  efficiency  of  the  pylons  is  shown  in  Fig.  13.  In  general  mixing  efficiency  increases  with 
downstream  distance.  The  alternating  wedge  pylon  exhibits  the  greatest  mixing  efficiency  of  all  the  pylons 
throughout  the  wake  region.  The  ramp  pylon  shows  a  reduced  mixing  efficiency  to  the  alternating  wedge 
pylon,  and  exceeds  the  mixing  efficiency  of  the  basic  pylon  in  most  of  the  wake  region.  The  basic  pylon  has 
a  reduced  mixing  efficiency  compared  to  both  the  alternating  wedge  and  ramp  pylons  until  the  very  end  of 
the  wake  region  computed. 

The  flammable  plume  extent  is  shown  in  Fig.  14.  The  alternating  wedge  pylon  has  the  largest  flammable 
plume  area  throughout  the  wake  region.  The  ramp  pylon  has  the  second  largest  flammable  plume  area 
in  most  of  the  wake  region,  however,  the  basic  pylon  flammable  plume  area  matches  the  ramp  pylon  as 
downstream  distance  increases.  The  flammable  plume  fractions  reveal  the  most  advantageous  location  to 
ignite  the  fuel  plumes.  The  alternating  wedge  and  ramp  pylons  both  show  a  rise  to  maximum  flammable 
plume  fraction,  10  de,  upstream  of  the  basic  pylon  rise  to  maximum  flammable  plume  fraction,  20  de. 
The  ramp  and  alternating  wedge  pylons  have  larger  fractions  of  the  fuel  plume  ready  for  combustion  in  a 
shorter  downstream  distance  than  the  basic  pylon.  Between  the  the  ramp  and  alternating  wedge  pylons,  the 
flammable  plume  fraction  of  the  alternating  wedge  rises  slightly  faster  in  the  near- field  region. 

A  last  comparison  is  to  examine  the  tradeoff  between  pressure  loss  created  by  each  pylon  versus  the 
potential  heat  energy  they  could  produce  through  mixing  the  injected  fuel.  This  comparison  is  made  at  each 
of  the  eight  positions  in  the  wake  region.  Equation  9,  potential  work  change,  is  the  method  of  comparison. 
At  each  axial  position  in  the  wake  region  the  mass  flow  of  ethylene  (rhf),  the  total  mass  flow  (m),  the  inverse 
total  pressure  loss  ratio  ((1  —  A)-1),  the  mixing  efficiency  (rym),  and  the  flammable  plume  fraction  (FPf) 
are  calculated.  These  measures  are  used  in  the  tradeoff  calculation.  The  other  variables  in  the  equation  are 
held  constant,  hf  =  4.5E7  J/kg:  this  is  a  typical  heating  value  for  a  hydrocarbon  fuel.  R  =  287  J/kg  K: 
although  ethylene  has  a  different  gas  constant,  it  is  in  the  minority,  so  the  gas  constant  of  air  is  used.  Tj 
=  712  K:  this  is  a  typical  combustor  inlet  static  temperature,  calculated  in  previous  work.1  The  ratio  of 
potential  work  change  between  the  ramp  and  alternating  wedge  pylons  compared  to  the  basic  pylon  is  shown 
in  Fig.  15. 
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The  alternating  wedge  and  ramp  pylons  both  exceed  the  basic  pylon  in  potential  work  change  for  a  large 
portion  of  the  wake  region.  It  is  noteworthy  that  the  basic  pylon  matches  the  ramp  and  alternating  wedge 
pylons  around  29  de.  At  this  point  the  two  hypermixer  pylons  loose  their  energy  advantage  over  the  basic 
pylon.  However,  the  flammable  plume  mixture  downstream  of  the  hypermixer  pylons  is  ready  for  combustion 
around  10  de,  and  they  still  have  the  energy  advantage  over  the  basic  pylon  at  that  location.  In  summary, 
there  is  a  finite  distance  downstream  where  the  alternating  wedge  and  ramp  pylons  perform  better  than  the 
basic  pylon.  Between  the  two  hypermixer  pylons,  the  alternating  wedge  pylon  out  performs  the  ramp  pylon 
slightly. 


VIII.  Comparison  to  Experiment 

A  qualitative  comparison  of  simulation  and  initial  experimental  data  obtained  is  presented.  There  are 
some  differences  in  conditions  between  simulation  and  experiment.  The  inviscid  walls  of  the  simulation 
versus  the  viscous  walls  of  the  experiment  is  one  difference.  Also,  the  primary  flow  Mach  number  is  2.0  in 
the  simulation  and  1.95  in  the  wind  tunnel  experiment.  The  ramp  and  alternating  wedge  pylon  simulations 
are  compared  to  experiment  since  they  have  a  reasonably  steady  wake  flow.  The  data  gathered  in  the  wind 
tunnel  was  done  using  a  static  pressure  cone  probe  and  a  total  pressure  pitot  probe,  yielding  averages  at  a 
measurement  point.  The  wind  tunnel  probing  location  is  6.44  inches  downstream  of  the  pylon  (23.9  de) ,  and 
the  probing  area  is  2  in  x  3  in  in  extent.  Figure  16  shows  the  probing  area  in  relation  to  the  pylon. 

Mach  number  is  the  parameter  chosen  for  comparison.  The  ramp  and  alternating  wedge  pylons  both  show 
similar  wake  region  shapes  and  Mach  number  magnitudes  in  experiment  and  simulation  (Fig.  17  and  18). 
The  resolution  of  the  simulation  data  is  much  better  than  experiment.  The  probing  area  holds  approximately 
3500  simulation  cells,  and  only  225  probing  locations.  One  notable  difference  in  the  Mach  profiles  is  the 
side  areas.  The  experimental  results  show  higher  Mach  numbers  on  the  sides.  This  is  most  likely  due  to  the 
reflection  of  the  shock  wave  off  the  wall  being  in  a  different  location  due  to  the  wall  boundary  layer  present 
in  the  wind  tunnel  experiment,  but  not  in  the  simulation.  The  position  of  the  reflected  shock  could  also 
account  for  the  slight  differences  in  the  shapes  of  the  pylon  wakes.  Overall,  the  Mach  profiles  are  similar  in 
simulation  and  experiment  for  these  two  pylons. 


Figure  16.  Wind  tunnel  probing  area. 
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(a)  CFD.  (b)  Wind  tunnel  data. 

Figure  17.  Ramp  pylon  Mach  number. 


120 

100 

80 

e" 

E_ 

>  60 

40 

20 


40  60  -60  -40  -20  0  20  40  60 

Z  (mm) 

(b)  Wind  tunnel  data. 

Figure  18.  Alternating  wedge  pylon  Mach  number. 
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IX.  Conclusions 

Three  pylons  are  compared  in  their  performance  as  in-stream  scramjet  fuel  injectors.  All  three  pylons 
exhibit  the  same  frontal  area  to  the  flow.  A  cold  flow  CFD  study  with  ethylene  fuel  injection  is  accomplished 
for  each  pylon.  The  ramp  and  alternating  wedge  pylons  show  decisive  near-field  wake  region  increases  in  axial 
vorticity  compared  to  the  basic  pylon.  They  also  exhibit  an  increase  in  total  pressure  loss  accompanying 
the  increased  streamwise  vorticity.  The  increased  vortical  motion  results  in  better  mixing  characteristics  for 
those  pylons,  bringing  the  ethylene/air  mixture  to  flammable  conditions  in  a  shorter  distance.  This  promotes 
the  idea  that  streamwise  vorticity  production  can  be  used  to  shorten  combustor  lengths.  A  theoretical  energy 
comparison  also  shows  that  the  ramp  and  alternating  wedge  pylons  have  an  energy  advantage  over  the  basic 
pylon  for  a  finite  distance  downstream  of  the  pylon  if  the  fuel  plume  is  ignited  within  this  distance.  A  last 
observation  is  the  alternating  wedge  pylon  has  a  slight  performance  advantage  over  the  ramp  pylon. 

The  views  expressed  in  this  article  are  those  of  the  authors  and  do  not  reflect  the  official  policy  or  posi¬ 
tion  of  the  United  States  Air  Force,  Department  of  Defense,  or  the  U.S.  Government. 
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